Thermoelectric properties of the bismuth telluride nanowires 
in the constant-relaxation-time approximation 



Igor Bejenarl3 

Institute of Electronic Engineering and Industrial Technologies, 
ASM, 3/3 Academiei str., MD2028 Kishinev, Moldova and 
Bogoliubov Laboratory of Theoretical Physics, 
Joint Institute for Nuclear Research, 
14 1980 Dubna, Moscow region, Russi^ 

Valeriu Kantsei0 
Institute of Electronic Engineering and Industrial Technologies, 
ASM, 3/3 Academiei str., MD2028 Kishinev, Moldova 
(Dated: August 18, 2008) 

Electronic structure of bismuth telluride nanowires with the growth directions [110] and [015] 
is studied in the framework of the anisotropic effective mass method using the parabolic band 
approximation. The components of the electron and hole effective mass tensors for six valleys are 
calculated for both growth directions. For a square nanowire, in the temperature range from 77 K 
to 500 K, the dependence of the Seebeck coefficient S, the thermal k and electrical conductivity g 
as well as the figure of merit ZT on the nanowire thickness and on the excess hole concentration pex 
are investigated in the constant-relaxation-time approximation. The carrier confinement is shown 
to play essential role for nanowires with cross section less than 30 x 30 nm^. In contrast to the 
excess holes (impurities), the confinement decreases both the carrier concentration and the thermal 
conductivity but increases the maximum value of the Seebeck coefficient. The confinement effect is 
stronger for the direction [015] than for the direction [110] due to the carrier mass difference for these 
directions. In the restricted temperature range, the size quantum limit is valid when the p— type 
nanowire cross section is smaller than 8 x 10 nm^ (6 x 7nm^ and 5 x 5 nm^) at the excess hole 
concentration p^x = 2 x lO^^cm"'^ {p^x = 5 x 10^*cm~^ and Pex = 1 x lO^^cm"^ correspondingly). 
The carrier confinement increases the maximum value of ZT and shifts it towards high temperatures. 
For the growth direction [110], the maximum value of the figure of merit for the p— type nanowire is 
equal to 1.4, 1.6, and 2.8, correspondingly, at temperatures 310 K, 390 K, and 480 K and the cross 
sections 30 x 30 nm^, 15 x 15 nm^, and 7x7 nm^ {pex = 5 x 10^* cm~^). At the room temperature, 
the figure of merit equals 1.2, 1.3, and 1.7, respectively. 

PACS numbers: 73.63.Nm; 73.50.Lw 



I. INTRODUCTION 



The study of nanowire systems is of interest be- 
cause of their possible applications, in particular, in 
high-efficiency thermoelectric devices. The strong two- 
dimensional confinement of such systems allows manip- 
ulating kinetic effects P, [3, The reduction of di- 
mensions leads: (i) to an increase in the Seebeck coef- 
ficient because of the increasing density of states in the 
vicinity of the Fermi energy, (ii) to a possible use of the 
anisotropy factor of the Fermi surface in the multi-valley 
semiconductors, (iii) to an increase in the phonon bound- 
ary scattering at interfaces of a heterostructure with- 
out an appreciable increase in electron boundary scat- 
tering, (iv) to an increase in the carrier mobility at a 
given concentration due to the size quantization effect. 
The most promising thermoelements for cooler manu- 
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facturing are those based on the nanostructures consist- 
ing of anisotropic materials like bismuth and lead tel- 
lurides [1, Q. Bismuth telluride is such a multi- valley 
material with a highly anisotropic isoenergy surface near 
the L— point in the Brillouin zone. The charge carrier 
effective mass anisotropy enhances the carrier confine- 
ment effect and complicates the carrier motion through 
a cross section of the nanowire. The anisotropy fac- 
tor leads to a modification of physical properties of the 
nanowire. For instance, the semimetal-semiconductor 
transition takes place at the Bi nanowire diameter de- 
pendent on the effective mass anisotropy [1, Bis- 
muth telluride and its solid solutions {Bi2-xShxTe^, 
Bi2Te3-ySey) are presently the best thermoelectric ma- 
terials for commercial applications at room temperature. 
A large value of the thermoelectric efficiency Z of these 
materials is due to the high degeneracy of the energy 
band edges. It is necessary to increase the Seebeck coef- 
ficient or the electric conductivity in order to increase 
the thermoelectric efficiency. An increase in the See- 
beck coefficient is provided by the increase in the carrier 
effective mass, while the electric conductivity increases 
due to a decrease in the effective mass. The effective 
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mass anisotropy factor allows for bypassing this problem 
and satisfying both conditions. Fast development of nan- 
otechnologies provided various techniques to prepare bis- 
muth telluride nanowires. Monocrystal and polycrystal 
BixTei-x nanowires with diameters 40-;-60nm are ob- 
tained by means of the electrochemical deposition in the 
nanopores of anodized alumina membranes Free- 
standing Bi and Bi2Tez nanowires are fabricated using 
both the Ulitovsky technique and the method of a high- 
pressure injection of the melt into capillaries 0. The di- 
ameter of such nanowires varies from 100 nm to 10 /im. 
There are only two growth directions [110] and [015] for 
bismuth telluride nanowires. 

The thermal conductivity k of Bi2Te^ nanowires with 
diameter 40 nm was experimentally shown to be reduced 
at least by an order of magnitude from the bulk value 
due to the phonon-boundary scattering, which dominates 
phonon-phonon Umklapp scattering at room tempera- 
ture ^3:]. The Seebeck coefficient S was measured to be 
by 15% -=-60% larger than the bulk values This 
thermal conductivity reduction and the aforementioned 
increase in the Seebeck coefficient in Bi2Te3 nanowires 
can be exploited for enhancing the thermoelectric figure 
of merit. There are only a few experimental results re- 
lated to thermoelectric parameters of bismuth telluride 
nanowires. 

Here, we consider thermoelectric parameters of the 
rectangular bismuth telluride nanowires with growth di- 
rections [015] and [110] using a constant-relaxation-time 
approximation, when both carrier energy and momen- 
tum are conserved in scattering processes. The carrier 
mobility coincides with the bulk value. The calculation 
method of kinetic coefficients in this case is similar to that 
treated by Lin for monopolar cylindrical Bi nanowires . 
We extended this method for an intrinsic semiconductor 
nanowire, taking into consideration the most of energy 
subbands in the bulk conduction and valence bands. The 
relaxation time depends on the carrier energy beyond the 
constant-relaxation-time approximation as r(e) = tqE^ . 
Fitting the experimental data to the theoretical expres- 
sion for the Seebeck coefficient gives the scattering factor 
r which varies from -1/2 to 1/2 in the bulk bismuth tel- 
luride material when 77K<T<300K [ll|- Therefore, 
our approximation, r — 0, is adequate for bismuth tel- 
luride nanowires when T>77K. 

Quantum well structures and nanowires based on 
Bi2Te3 have been studied previously in the size quan- 
tum limit (SQL) taking into account only the lowest sub- 
band using the bulk-effective-mass approximation 0, IT^ ] . 
The nanowire growth direction was supposed to be along 
one of the crystallographic axes. We recall that the effi- 
ciency of a thermoelectric material is characterized by 
means of the figure of merit ZT — aTS^/K, where a 
is the electrical conductivity. The maximum value of 
the figure of merit was obtained to be equal to ZT = 5 
for a 5A-thick quantum well and ZT = 14 for a 5A- 
wide quantum wire. ZT significantly increases when the 
square nanowire width becomes smaller, by the order 



of magnitude, than the thermal de Broglie wavelength 
Xd « {h?/2mkBTf/'^. We wih show that the SQL is 
valid only for nanowires with very small cross section 
(8 X 8nm^) in a restricted temperature range. To treat 
the nanowires with larger cross sections at temperature 
77 K and higher, we consider all subbands of the conduc- 
tion and valence bands. For the calculation, we mainly 
use conventional bulk Bi2Te3 parameters with the ex- 
ception of the lattice thermal conductivity, which is less 
by an order of magnitude than the bulk value. Never- 
theless, we suppose the temperature dependence of the 
phonon thermal conductivity to be similar to that for 
the bulk material because of the absence of correspond- 
ing experimental data. The carrier confinement leads 
to a splitting of electronic band structure into subbands, 
while the phonon dispersion is assumed to be unchanged. 
The boundary scattering is not taken into account in our 
calculations. In general, we have improved the values 
of thermoelectric parameters of Bi2Te3 nanowires ob- 
tained with the constant-relaxation-time approximation 
and provided a background for more sophisticated calcu- 
lations with account of different scattering mechanisms. 

In the next section, we consider the longitudinal and 
cross-sectional effective mass components as well as the 
electronic band structure of six carrier pockets in the 
parabolic-band approximation for the rectangular bis- 
muth telluride nanowires with growth directions [015] 
and [110]. The parabolic-band approximation is ade- 
quate for a calculation of the electron energy spectrum 
because bismuth telluride is an indirect semiconductor. 
We take into account the temperature dependence of the 
effective mass of the carriers, which is appropriate for 
the intrinsic bismuth telluride. The temperature depen- 
dence of the electron and hole subbands is also consid- 
ered to obtain more realistic values of the temperature 
dependent thermoelectric parameters. We analyze condi- 
tions at which the carrier confinement effect on transport 
properties of bismuth telluride nanowires can be experi- 
mentally fixed. We also consider the nanowire cross sec- 
tion and the excess hole concentration at which the size 
quantum limit is adequate. In Sec. 3, a dependence of the 
Fermi energy and the carrier concentration on tempera- 
ture is studied for the nanowires with different cross sec- 
tions and the excess hole concentrations. In Sec. 4, the 
temperature dependence of the Seebeck coefficient, the 
thermal conductivity, and the figure of merit are stud- 
ied. Finally, conclusions are given in Sec. 5. 



II. ELECTRONIC STRUCTURE 
A. Electronic Band Structure of the Bulk Material 

The crystal structure of bismuth telluride with rhom- 
bohedral unit cell belongs to the symmetry group D^^^ 
(mm) [H, m. The crystallographic axes are binary 
n (x), bisectrix s (j/), and trigonal c (z). The electronic 
band structure of the indirect band gap Bi2Te3 semi- 
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conductor is given by a parabolic-band approximation 
in the Drable-Wolf six-valley model. All six valleys are 
equivalent in the bulk. The tilt angle between the hole 
(electron) energy ellipsoid principle axes and bisectrix- 
binary plane is 32° (34°) [l5|, [iGj.The measurement of 
the Shubnikov-de Haas effect showed the bulk effec- 
tive mass tensor components to be to/ji = 0.0308 mo, 
= 0.442 mo, and to^3 = 0.0862 mo for holes and 
niei = 0.0213 mo, me2 — 0.319 mo, and me3 = 0.0813 mo 
for electrons in the local coordinate system of the ellip- 
soid. The energy band gap was obtained by means of 
the measurement of the absorption coefficient and the 
electrical conductivity. The indirect (direct) band gap 
is equal to 0.15 eV (0.22 eV) at temperature 2 K. It is 
0.13 eV at room temperature, being about linearly de- 



pendent on temperature with the rate dEg/dT 



-0.09 



meV/T. This agrees rather well with the recent calcu- 
lations made within the screened-exchange local den- 
sity approximation [l7[. Using the general formula 
Eg{T) = Eg{fi) ~ aT'^/iT + (3) reliable for most semicon- 
ductors and the experimental data listed above we ob- 
tained the following expression for temperature depen- 
dence of Bi2Tez energy band gap [H 



£;,(T) = 150- 



0.0947 r2 
T-f 122.5 



[meV]. 



(1) 



The electrical conductivity measured in the bis- 
muth telluride bulk material varies with temperature 
as (j(T) « 7^-2 _^ rp-i.2 Tffi^QTi^ carrier concentration 
is about lO^^cm-^ [l|, [l^, [SI. For classical statis- 
tics, when the acoustic phonon scattering is dominant, 
the electrical conductivity depends on temperature as 
cr(T) « ^] Hence, in the case of the in- 

trinsic conductivity, the temperature dependence of the 
effective mass is given by mh ~ T^-^'^ T^-^ for holes and 
me « T°-^^ for electrons [23| . 

The parabolic valence band structure is described by 
the Schrodinger equation in the rhombohedral coordinate 
system for the bulk bismuth telluride as follows [T3| 



y VaV^'(r) = £;*(r), 



(2) 



where a is the inverse tensor of the hole effective 



mass with components an — 32.5 mp ^, a22 = 4.81 : 



a33 = 9.02 mn 



and a23 = 4.15 mg 



[15|. 



The 



edge of the Bi2Te^ valence band is displaced from 
the edge of the conduction band by the vector 
ko = (0.091, 0.152, 0.152) nm^^ in the momentum space 
[l3]- Due to this displacement, the conduction band 
structure is described as 



y(-^V-ko)a(-^V-ko)*(r) 



E-^{v). (3) 



The components of the inverse electron effective 



mass 



tensor are aii=46.9mg^, a22 = 5.92 mg 



a33 = 9.50 mg ^, and a23 = 4.22m^i [3. Scattering of 
charge carriers from different valleys {X L) is simi- 
lar to that from equivalent valleys, the difference of the 



carrier effective masses and the edges of different valleys 
being taken into account Hence, the displacement 

vector ko of the edges can be neglected in the calculation 
of kinetic coefficients. 



B. Anisotropic Band Structure in the Bismuth 
Telluride Nanowire 



The growth direction [110] of a Bi2Te3 nanowire with 
diameter 50 100 nm is perpendicular to the trigonal 
axis. This direction is equivalent to [1120] in the hexag- 
onal unit cell. It is parallel to one of three bisectrix axes 
perpendicular to the trigonal axis Q. We choose the 
nanowire coordinate system with y axis along the bisec- 
trix axis to apply the boundary condition to Eqs. ([2]) and 
We suppose that the wave function vanishes at the 
nanowire boundary. Two direct equivalent valleys 1 and 
4 are distinguished from the set of four oblique equiv- 
alent valleys 2, 3, 5, and 6 for the [110] growth direc- 
tion. Since the carrier motion in the y-direction is free, 
^1,4(1") = u(x, z) exp[z/cj;(j/ — za23/a33)] is an appropri- 
ate wave function of the hole pockets 1 and 4 aligned 
with the bisectrix axis and Eq. ^ can be rewritten as 
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an 



"33 



92 
dz 



= E 



b 

2ml 



u{x,z), (4) 



where the longitudinal effective mass component is 



(a22 



^^23 



/a33) 



The transverse effective mass 



components are m^, = l/an and mz = 1/033. 

For a self-consistent calculation of the energy spectra 
for different valleys, we use a (rhombohedral) coordinate 
system associated with the nanowire. Therefore, one has 
to rotate the hole (electron) pockets 2, 3, 5, and 6 about 
z (c) axis to obtain the corresponding inverse effective 



mass tensor d 



3(P) 



R 



in the rhombohedral 



system. The rotation operator is 



cos[(p- l)7r/3] sin[(p - l)7r/3] 
— sin[(p — l)7r/3] cos[(p — l)7r/3] 
OA 



(5) 



The parameter p indicates the number of the correspond- 
ing valley (1, 2,..., 6); A = +l ffp = 1,2,3 and A = -1 
li p ~ A, 5, 6. To eliminate the mixed partial derivative 
from the Schrodinger equation for the oblique valleys, 
one has also to transform the rhombohedral coordinate 
system into a new one {x',y',z') by using the rotation 
operator R^z about y (or s) axis by the angle 8 defined 
from the relation tan(2e) = ±2a-^^ / {a^^ - a^^^). In the 
new coordinate system, the wave function for the oblique 
valleys is 



F(P)(r') 



u[x , z ) eyi'pyikyy) exp [ —i — ^^^"^ 



X exp \ -i^^^z' 



(6) 
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A function u{x' , z') satisfies the equation in form 
(PI) with a transformed inverse effective mass tensor 
p{p) — R~la^P^ Rxz- Therefore, the longitudinal effective 
mass component for the oblique valleys 2, 3, 5, and 6 is 



(7) 



The transverse 



effective 

(p) 



mass components are 



rnf^ = f//3j^' and to 

BioAeTeo.bi nanowires with diameter 40 60 nm 
mainly grow along the [Of 5] direction This direction 
lies along the vector 1 = (0, 0.8489, 0.5285) in the rhom- 
bohedral coordinate system. In general, the growth di- 
rection of the bismuth telluride nanowire tends towards 
the trigonal axis with decreasing diameter in the same 
way as a Bi nanowire does 0. To satisfy the bound- 
ary conditions, Eqs. ([2]) and Q should be written in a 
nanowire coordinate system (x" , y" , z" ) with the y" axis 
along the vector 1. We have to use the rotation matrix 
Ryz to transform the rhombohedral coordinate system 
into the nanowire coordinate system. The rotation is 
about the binary axis by the angle 8' = arccos(0.8489) 
(32°). Hence, this matrix has the form 



Ryz — 



1 
0.8489 
0.5285 





-0.5285 
0.8489 



(8) 



In the new system of coordinates, the wave function and 
the longitudinal effective mass component for six valleys 
are written in the form ^ and (O, correspondingly. The 
components of the inverse effective mass tensor are de- 
fined as 7^^' = R~}R~}a'^P^ RyzRxz- In this case, the ro- 
tation angle Q is slightly modified taking into account 
the transformation Ryz- 

The transverse and longitudinal hole (electron) effec- 
tive mass components in the bismuth telluride nanowires 
oriented along the bisectrix axis and the direction [015] 
are listed in Table U ([IT]). The hole effective mass is 
greater than that of electron both in bulk and in the 
nanowire. For the nanowire growth direction [015], the 
carrier effective mass components for the 1st and 4th val- 
ley slightly differ from the bulk components because this 
direction is approximately oriented along one of the prin- 
ciple axis of the bulk electron (hole) ellipsoid. Since the 
effective mass components are alike for valleys 2, 3, 5, and 
6 in the nanowire oriented along the bisectrix axis, these 
valleys are equivalent in this case. For the nanowires 
with the growth direction [015], some of the hole (elec- 
tron) effective mass components in valleys 2 and 5 differ 
from the corresponding components in valleys 3 and 6 
approximately by a factor of four (eight). The transverse 
effective mass anisotropy of carriers in valleys 1 and 4 
(3 and 6) for the nanowires oriented along the bisectrix 
axis is greater (less) than that for the nanowires with the 
growth direction [015]. The mass anisotropy for valleys 
2 and 5 is similar for both growth directions. For valleys 



1, 2, 4, and 5, the cyclotron effective mass ^JraxTfiz for 
the direction [015] is less than that for nanowires oriented 
along the bisectrix axis. For nanowires with the growth 
direction [015], the effective mass anisotropy in valleys 3 
and 6 is greater than that for the other direction. Due to 
last two factors, the confinement effect for the nanowire 
with the growth direction [015] is larger than that for 
the nanowire with the growth direction along the bisec- 
trix axis. 

Here, we study bismuth telluride nanowires with a rect- 
angular cross section. Similar bismuth nanowires have 
been recently fabricated using the electron beam lithog- 
raphy |2^ . The energy spectrum of holes is defined from 
Eq. (g]) as follows 



rrhole 



(ky 



ruxal mzol 



2TOr, ' 



(9) 



The electron energy is defined from Eq. ^ in the form 



2 



{n-K/ax-koxf , {lir/az - kozf 



h^{ky - kpy) 
2m*, 



(10) 



Symbols and Uz denote the sides of the rectangular 
cross section. The carrier effective mass components rux , 
TO*, and niz are listed in Tables HI and HIl The electron 
and hole wave functions differ by the oscillating factor 
exp(zkor) which does not effect the transport matrix el- 
ements. 
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FIG. 1: Electron and hole energy spectrum for valleys 1 and 
4 in the intrinsic nanowire {a^ = ttz = 15nm) with growth di- 
rection [110]. The solid lines correspond to quantum numbers 
n = 1,2 and I — 1, dotted lines to Z = 2, dashed lines to Z = 3. 

Fig. [T] shows the electron and hole energy spectrum in 
the intrinsic nanowire with the growth direction along 
the bisectrix axis at room temperature. The spectrum is 
composed of six subbands for equivalent valleys 1 and 4. 
The sides of the nanowire cross section are a^; = az = 15 
nm. The energy origin is taken at the top of the va- 
lence band in the bulk. The bottom of the electron sub- 
bands is displaced from the top of the valence subbands 
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TABLE I: Calculated effective mass components of a hole pocket in the bismuth telluride nanowire along the indicated crys- 
tallographic directions. The y direction is chosen along the nanowire growth direction (the wire axis). All values of the mass 
components are in units of the free electron mass mo. 







Bisectrix 






[015] 




Hole mass 


1st and 4th 


2nd and 5th 


3rd and 6th 


1st and 4th 


2nd and 5th 


3rd and 6th 


component 


pocket 


pocket 


pocket 


pocket 


pocket 


pocket 




0.03077 


0.07034 


0.07034 


0.03077 


0.1017 


0.4220 


ml 


0.3448 


0.1093 


0.1093 


0.4421 


0.2013 


0.05893 


mz 


0.1109 


0.1530 


0.1530 


0.08645 


0.05744 


0.04729 



TABLE IL Calculated effective mass components of an electron pocket in the bismuth telluride nanowire along the indicated 
crystallographic directions. The y direction is chosen along the nanowire growth direction (the wire axis). All values of the 
mass components are in units of the free electron mass mo. 







Bisectrix 






[015] 




Electron mass 


1st and 4th 


2nd and 5th 


3rd and 6th 


1st and 4th 


2nd and 5th 


3rd and 6th 


component 


pocket 


pocket 


pocket 


pocket 


pocket 


pocket 


rrix 


0.02132 


0.05625 


0.05625 


0.02132 


0.08977 


0.03550 


ml 


0.2472 


0.07779 


0.07779 


0.3197 


0.1484 


0.04982 


mz 


0.1053 


0.1268 


0.1268 


0.08139 


0.04166 


0.3137 



along the axis ky by 0.152 nm. Fig. 1 demonstrates that 
the size quantization is manifested more significantly for 
holes than for electrons. However, such a large difference 
between the electron and hole confinement is an apparent 
effect. The real difference is smaller. The sketch in Fig. 1 
is reflected when the energy origin is situated at the bot- 
tom of the bulk conduction band, while the confinement 
effect does not depend on the choice of the energy origin. 
A slight difference between the electron and hole con- 
finement exists because the electron and hole masses are 
not equal. The energy separation between the subbands 
with different quantum numbers n is greater than that 
between the subbands with different quantum numbers I 
due to the ratio of the transversal mass components rux 
and ruz- The numerical calculations demonstrate that a 
proper choice of the sides of the rectangular cross section 
allows for manipulating the energy subband splitting. A 
fitting of the ratios of the sides diminishes or enhances 
the mass anisotropy effect on the energy spectrum. This 
agrees with the conclusion of our previous research on Bi 
nanowires We have thus obtained that the model of 
carriers with an anisotropic mass in the nanowire with a 
circular cross section is identical with the model of carri- 
ers with an isotropic mass in the nanowire with an elliptic 
cross section. 



Evaluation Criteria of a Quantum Size Effect in 
the Bismuth Telluride Nanowire 



An experimental observation of the quantum size effect 
is possible at some favorable conditions. These conditions 
are mainly based on three criteria. The first criterion 
states that a thermal excitation should be small enough 
to provide interband carrier transitions, which diminish 



the effect of the band splitting on the transport prop- 
erties. In the so-called "classical limit", if the thermal 
energy fc^T is much larger than the energy level spac- 
ing AE, the subbands can be treated as continua owing 
to the thermal smearing of states 26[. In the "quan- 
tum regime" , variations in the current peak height versus 
gate voltage measured, for example, in Si nanowire de- 
vice suggest that the transport is carried out through co- 
herent energy states with the energy level spacing larger 
than the thermal energy p7| . The subband splitting in- 
creases with increase of quantum numbers in the model 
of a cylindrical potential well with a rectangular cross 
section. Hence, taking into account Eq. ^ , the first cri- 
terion can be mathematically expressed as 



max[E2.i{ky), Ei,2iky)] - Ei,i{ky) > fc^T. 



(11) 



We choose the maximum energy values of subbands with 
quantum numbers (1,2) and (2,1) in the relation (fTTj) 
because the transition matrix elements including the re- 
spective wave functions Fi2(r) and F2i(r) are slightly 
different. In the case of a degenerate semiconductor, the 
quantum size effect is reduced when there are many elec- 
tron subbands under the Fermi level. Since the carrier 
thermal excitation is much less than the Fermi energy, 
the first criterion for a degenerate nanowire reads 



End + l{ky) — En,l{ky) > Ep, 



(12) 



where the quantum numbers n and I are small. Condi- 
tions pT|) and (fT2l) are equivalent to the requirement that 
the de Broghe wavelength {Xd) of electrons should be of 
order of or larger than the nanowire thickness a [2^ 



A, 



> 1. 



(13) 
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In the case of classical statistics, the de Broglie wave- 
length is defined as Xu ~ 2nh/ ^/2mkBT and the condi- 
tion (jlip can be obtained from Eq. (fT3|) . Taking into 
account that Xd ~ Xf ^ 2'Kh/ \/2mEp for a degenerate 
semiconductor nanowire, the condition (|13[) can be trans- 
formed into the expression JT/ 




FIG. 2: Temperature dependence of the reduced energy sepa- 
ration between the hole subbands (1,1) and (2,1) in the square 
nanowire, with thickness a = 30 nm (dashed-dotted line), 15 
nm (dashed line), and 7 nm (solid line) for the 1st and 4th 
valley. 

Fig. [2] demonstrates the temperature dependence of 
the energy separation between the hole subbands with 
quantum numbers (1,1) and (2,1) in the units of ther- 
mal energy in the square nanowires with growth orien- 
tation [110]. The subband separation in the nanowires 
with thickness 30 nm tends to unity when T > 350 K. 
The first criterion is satisfied in the nanowires with 
thickness less than 30 nm for whole temperature range, 
77 K < r < 500 K. Hence, the carrier confinement af- 
fects the transport properties in the nanowires with cross 
section less than 30 x 30nm^. The results for both 
nanowire growth directions are approximately similar. 

An electron (hole) can occupy the quantum state dur- 
ing a certain interval of time t because of scattering. The 
final lifetime of the quantum state leads to uncertainty 
of the corresponding energy due to the energy-time un- 
certainty principle. Energies E and E + AE cannot be 
resolved if the mean time between collisions r is shorter 
than that given by the energy-time uncertainty AEt > h 
because of a collision broadening of the energy levels by 
AE [2^. Therefore, the second criterion is that the space 
between the discrete energy levels must be much greater 
than the energy uncertainty, in order to keep the peculiar- 
ity of the one-dimensional density of states in nanowire. 
It is mathematically expressed for a rectangular nanowire 
as 



E2,i — El l 3> 



(14) 



h he 

To estimate the criterion, we use the value of 
the electron (hole) mobility /ig^^"' = 1200 cm^/Vs 
(/^L^^''^ =510 cm^/Vs) obtained from the measurement 



of the electrical conductivity in the cleavage plane of a 
5^2^63 monocrystal at 300 K [H, [l3, [ll . The temper- 
ature dependence of the electron (hole) mobility is given 
by the relation /ie « T"! "^ (/i/^ « T-^-O) [ll|, M, 
Using the mass components listed above, we estimate 
the electron and hole density-of-state effective mass 

TTiDOS = -^Jeg (™i™2TO3)^^'^ for the six-fold degenerate 
{Ndeg — 6) bulk valleys to be equal to 0.271 toq and 0.348 
mo, respectively. Hence, the electron (hole) energy un- 
certainty h/r varies from 6.5 meV (3.6 meV) to 16 meV 
(10 meV) in the temperature range 300 K <T < 500 K. 
The temperature dependence of the effective mass leads 
to a shrinkage of limits of the energy uncertainty inter- 
val. The electron (hole) thermal excitation ranges from 
26 meV to 43 meV in the same temperature interval. 
Therefore, the second criterion ([H]) is always satisfied 
to observe experimentally the confinement effect in the 
bismuth telluride nanowires, when the first criterion (jlip 
or in]) is fulfilled. 

In the ballistic region, the electron scattering at the 
surface roughness of a quantum wire destroys the con- 
ductance quantization |32] . The third criterion requires 
the magnitude of roughness of the nanowire boundary to 
be less than the de Broglie wavelength of electrons, in or- 
der to diminish the role of the surface diffusion .sj] . 
The wavelength lower limit is Xd ==21.4 nm (13.6 nm) for 
the electrons and Ad = 19.5 nm (11.7 nm) for the holes 
at 300 K when the effective mass components are 0.1268 
mo (0.3137 mo) and 0.153 mo (0.422 mo), respectively, 
for the growth direction [110] ([015]). This result agrees 
well with the recently estimated thermal de Broglie wave- 
length of holes (Ad = 11.4 nm) in the Bi2Te3/ Sb2Te3 su- 
per lattices [15. Since the surface roughness of Bi2Te3 
nanowires (films) is of the order of 1 nm (10 nm), the 
third criterion is satisfied [1, [s^ . 

The three criteria considered above include both the 
carrier effective mass and the nanowire cross section area. 
To study the size quantum limit, we also have to take into 
account the excess hole concentration. SQL is achieved 
when all electrons (holes) occupy the lowest (highest) 
subband in a degenerate semiconductor nanowire. This 
condition is mathematically expressed, starting from the 
relation (fT2l). as 



niD < 



8m*JE2i - En) 



nh 



(15) 



where the right hand side represents a full number of 
states in the 1st subband per unit length, Gid{E2i)- For 
the p— type degenerate bismuth telluride nanowires, the 
hole concentration coincides with the excess hole con- 
centration Pel- Using the relation (|15p . we obtain the 
maximal nanowire cross section thickness and width 
Uz at which the SQL takes place, from the following re- 
lation 



12m; 

2 



1/6 



(16) 
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For nanowires with the growth direction [110], the max- 
imum thickness and width is 8 nm (7 nm) and 10 nm 
(8 nm) at a low temperature when the excess hole con- 
centration is 5 X 10^^ cm^^ (1 X 10"^^ cm~'^). For the 
nanowire growth direction [015], we estimate the max- 
imum thickness and width as 6 nm (5 nm) and 9 nm 
(7 nm), correspondingly. Hence, the SQL for nanowires 
grown along [110] with a small excess hole concentration 
is larger than that for the growth direction [015] because 
the longitudinal effective mass component of holes m* 
in the highest valence subband for the direction [110] is 
twice as large as that for the other direction (for valleys 3 
and 6). At a high temperature, the maximal cross section 
sides depend on temperature when the hole concentration 
deviates from the excess hole concentration. In the next 
sections, we study thermoelectric properties of nanowires 
with the cross sections 7x7 nm^ and 15x15 nm^ to con- 
sider both the size quantization effect and the size quan- 
tum limit for temperatures ranging from 77 to 500 K. 



III. CARRIER STATISTICS IN THE 
NANOWIRE 

Carrier concentration for an electron pocket in the 
nanowire is defined as [13, 



E 

1=1 



*-l/2(?7«,i)- 



(17) 



,1/2 



The factor N^^ = {iml^^kBT /irh^Y' ^ denotes the effec- 
tive density of states for the corresponding electron (hole) 
pocket in the nanowire. 



1 



T{j + l)J exp(e - 7?) -f 1 



dt 



(18) 



is the complete Fermi-Dirac integral with a fractional in- 
dex. The reduced chemical potentials 7] for an electron 
subband and a hole subband are defined as 



Ey + En,i ~ 



(19) 
(20) 



where ■* determines the electron (hole) subband edge 
in the nanowire relative to the bottom (top) of the con- 
duction (valence) band for the bismuth telluride bulk ma- 
terial. 

The Fermi energy Ep is calculated from the equation of 
the electrical neutrality n{Ep) = p[Ep) — pex, where Pex 
denotes the concentration of uncompensated acceptors 
(excess holes). The excess holes are absent in an intrin- 
sic nanowire. The calculation of the bismuth telluride 
conduction (valence) band structure in the screened- 
exchange local density approximation demonstrated that 



its width Ebw varies in the range 0.8 0.9 eV [l7!l. We 
suppose that all electron (hole) subbands in the nanowire 
are fitted into the conduction (valence) band with width 
corresponding to the bulk material. In fact, we can take 
Ebw — 0.6 cV because it is greater than 10 fc^T. The 
maximum values of the subband quantum numbers Umax 
and Imax included in Eq. (|17p are defined from the fol- 
lowing expressions 



^max (^) 



2Ei,wmxal 
2Eb 



rux 


( 0,x 


)] 


1/2 














' m 








mx J 



1/2 



(21) 



(22) 



The bulk BixTei^x material is of p— type for x > 0.37 
and of n— type for x < 0.37 ^sj. The dependence of the 
nanowire conductivity type on the Bi-to-Te ratio differs 
from that of the bulk [3,]. For the sake of simplicity, 
we suppose that the -Bjo.37Teo.63 nanowire is intrinsic. 
This does not actually influence the results because the 
parameter x is not taken into account in the calculation. 

Fig. [3] represents a dependence of the Fermi energy 
Ep and the carrier concentration Ui = nip)/ (axCiz) on 
temperature in square intrinsic -Bio. 37^60. 63 nanowires 
with thickness 15 nm, 7 nm and the growth directions 
[110] and [015]. The energy origin is selected in the mid- 
dle of the band gap. The carrier confinement increases 
the Fermi energy in the nanowire because the hole mass 
is greater than the electron mass. The increase of the 
Fermi energy is greater for the nanowire growth direc- 
tion [015] compared with the other direction at the same 
nanowire thickness. The difference between the Fermi 
energies associated with the directions [015] and [110] 
increases when both the nanowire cross section and tem- 
perature decrease. For the -810.377^60.63 nanowires, the 
Fermi energy lies in the band gap in the whole tempera- 
ture range. Therefore, the carrier confinement does not 
change the conductivity type of such nanowires. Hence, 
we take into account the temperature dependence of the 
effective electron (hole) mass mg = ('^^ _ 2^0.17^ 

in our calculation of transport properties of the intrin- 
sic bismuth telluride nanowires. This mass temperature 
dependence leads to an additional increase in the carrier 
concentration with temperature. Fig. 3(b) demonstrates 
that the confinement decreases the carrier concentration 
in the intrinsic nanowires. This agrees with the calcula- 
tions for the Bi nanowires . The confinement is greater 
for the direction [015] than for [110], therefore, the car- 
rier concentration for the direction [015] is somewhat less 
than for the other direction at the same value of the 
nanowire cross section. This small difference increases 
when temperature increases or the nanowire cross sec- 
tion decreases. 

Fig.|3]demonstrates the temperature dependence of the 
Fermi energy and the hole concentration for the p— type 
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T(K) T(K) 



FIG. 3; The temperature dependence of (a) the Fermi energy and (b) the carrier concentration in the square intrinsic 
Bio. 37^60. 63 nanowire with thickness 7 nm (sohd line), 15 nm (dashed hne) and the growth direction [110]. Dashed-dotted and 
dotted hnes correspond to thickness 15 nm and 7 nm for the growth direction [015]. FiUed squares show the Fermi energy and 
the carrier concentration in the bulk. 




FIG. 4: Temperature dependence of (a) the highest hole subband (filled squares) and the Fermi energy for the square type 
BiiTes nanowire with thickness 15 nm (dashed line) and 7 nm (solid line) at pex — 5 x 10^* cm~^ and pex = 1 x 10^^ cm~^ 
(dotted line) for the growth direction [110], (b) the hole concentration in the bulk (filled squares) and in the nanowire with 
thickness 7 nm (solid line) and 15 nm (dashed line) for the growth direction [110] at pex = 5 x 10^* cm~'^. Dashed-dotted 
(dotted) hne shows the Fermi energy (hole concentration) in the 7 nm (15 nm) thick nanowire with the growth direction [015] 
at Pex ^5 X 10^^ cm"^. 



Bi2Te3 nanowires with the cross section 15 x 15 nm^ 
and 7 x 7nm^ for the both growth directions. For the 
nanowire with the cross section 7x7 nm^ grown along 
the [110] direction, the highest hole subband of the 
equivalent valleys 2 and 6 is depicted by filled squares 
in Fig. mja). The energy origin is taken at the mid- 
dle of the band gap. The bulk Fermi energy varies 
from -100 mcV to -37 meV in the temperature inter- 
val 77 K < T < 500 K at the excess hole concentration 
Pex = 5 X 10^^ cm~^. If the nanowire cross section area 
decreases by a factor of five, the Fermi energy decreases 
by 40 per cent (50 per cent) at low (high) temperature. 
For the square nanowire with thickness 7 nm and the ex- 
cess hole concentration pex — 5 x 10^* cm~'^, the Fermi 
energy corresponding to the [110] and [015] growth di- 
rections lies under the highest valence subband when 
T < 180 K and T < 400 K, respectively. For the excess 
hole concentration pex = 1 x 10^^ cm~'^, the Fermi en- 
ergy is under the highest hole subband in the whole tem- 



perature range. The hole concentration in the nanowire 
depends neither on temperature nor on the cross sec- 
tion area at a high value of the excess hole concen- 
tration = 1 X 10^^ cm~'^. For a smaller excess hole 
concentration p^x = 5 x 10^* cm^'^ the hole concentra- 
tion in the nanowire with the cross section 15 x 15 nm^ 
slightly increases with temperature while the hole con- 
centration for the cross section 7x7 nm^ is constant 
in the whole temperature range. Therefore, the size 
quantum limit may be applied for consideration of the 
nanowires with the cross section 7x7 nm'^ at the ex- 
cess hole concentration pex = 1 x 10^^ cm^'^ in the whole 
temperature range, while it is valid only at low tempera- 
tures for a smaller excess hole concentration, for instance, 
Pex — 5 X 10^^ cm~^. 

The Fermi energy for the growth direction [015] is 
less than that for the [110] direction, nevertheless, the 
hole concentrations for the both nanowire growth direc- 
tions coincide at the cross section 7 x 7nm^. In con- 
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trast, the hole concentration for the growth direction 
[015] is somewhat less than that for the other direc- 
tion at pcx — 5 X 10^* cm"^, while the Fermi energy is 
approximately the same for the nanowire cross section 
15 X 15 nm^. The difference in the Fermi energies and the 
carrier concentrations associated with different growth 
directions is owing to the greater confinement effect along 
the [015] direction. Ep for the p— type nanowire signif- 
icantly changes with temperature compared to that for 
the intrinsic nanowire. The p— type nanowire hole con- 
centration slightly depends on temperature in compari- 
son with the intrinsic nanowire hole concentration. Both 
the carrier confinement and the excess holes decrease the 
temperature dependence of the nanowire carrier concen- 
tration; as a result, the intrinsic type of conductivity is 
suppressed in the p— type Bi2Te3 nanowires at high tem- 
peratures. Hence, the temperature dependence of the 
hole (electron) mass can be neglected in the transport cal- 
culations. In general, the size quantization decreases the 
carrier concentration for any type of conductivity of the 
nanowire. Our calculations show that the confinement 
effect increases the difference between the Fermi energies 
as well as between the carrier concentrations correspond- 
ing to the growth directions [110] and [015], in contrast 
to the effect of the excess holes. This leads to a differ- 
ence in the transport properties of the bismuth telluride 
nanowires with different growth orientations considered 
in the next section. 



IV. THERMOELECTRIC-RELATED 
TRANSPORT COEFFICIENTS 



A. Constant-relcLxation-time approximation 



The carrier mobility fi and the relaxation time 
T = To = nifi/e are supposed to be constant in the 
constant-relaxation-time approximation. The electrical 
conductivity a, the Seebeck coefficient S, and the elec- 
tron (hole) thermal conductivity Ke(ft) for a nanowire are 
defined as 0, [sl] 



s ^ - 



1 £(1) 

eT m ' 



e{h) 



1 



L(2)_ 



(23) 
(24) 

(25) 



We use the generalized transport matrix element L^") 
that presents the sum L^") — X)^!"'' over all six valleys. 

i 

The transport matrix elements for the ith valley have the 

form [rl-lsa] 



(0) 



n=l 1=1 

"max 'max(n) r 

LfUifcsTAE E i*i/2(<,)-<i'i>-i/2«i) 

n=l ;=1 L 



if) - (kBT^D, E E 

rt=l 1=1 



(26) 
(27) 
(28) 



where the sign " + " (" — ")corresponds to electrons 
(holes). The symbol Di denotes the maximal conduc- 
tivity of a nondegenerate nanowire associated with the 
ith valley, which is expressed as Di — eiiyN^^^^-^/axdz- 
The electron (hole) thermal conductivity can be written 
in the form Kf,(h) = '^e(h)LT , where the Lorenz number L 
for both degenerate one-dimensional and bulk semicon- 
ductor compounds is given by L'^^f = (7rfcB/e)^/3 [soj . 



For a bipolar nanowire, the transport matrix element 
includes both electron and hole parts, L^") = L" + L'^. 
The bipolar Seebeck coefficient Sj, and the thermal con- 
ductivity K are expressed as [H, 



S,^ + (29) 

K= Kl + Ke + + Keh, (30) 

where the last term k^h = Ta^crhiSh — SeY /{<Je + (^h) is 
attributed to the thermal energy carried by the electron- 
hole pairs generated at the heated end of the sample and 
moving to the cold end, where they annihilate. The gen- 
eration and annihilation of the electron-hole pairs are as- 
sociated with the absorption and emission of the thermal 
energy, correspondingly p^ . 

Eq. (29) and ((30)) can be reduced to those for monopo- 
lar semiconductor at the appropriate conditions. The 
first term kl in the right hand side of Eq. ([50]) is re- 
lated to the lattice thermal conductivity. At the room 
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temperature, the lattice thermal conductivity measured 

and that in the perpendic- 



along the trigonal axis k^"^' 



ular direction are equal to 0.725 W/mK and 1.45 

W/mK, respectively, for the bulk bismuth telluride mate- 
rial (4ll.l4^ . It is inversely proportional to temperature, 
KL = 430/r, when T > 50 K 43]. For the [015] direction, 
the lattice thermal conductivity is given by the expres- 
sion K^"^^' = cos^(6'')k^"^'^' -f- sin^(6'')«;|"'"^', where 9' de- 
notes the ang le between the bisectrix axis and the direc- 
tion [015] [20]. This conductivity is equal to 1.25 W/mK 
at the room temperature. Both theoretical estimations 
and the measurements of the lattice thermal conductiv- 
ity of the Bi2Te3 nanowire with diameter 40 nm indicate 
that the nanowire lattice thermal conductivity is reduced 
by an order of magnitude from the bulk value [l^, . 
We take into account this fact in our calculations. 

While the anisotropy factor of the thermal conductiv- 
ity is 2, the anisotropy factor of the electrical conductiv- 
ity is 4 (2.7) for both the n(p)— type bismuth telluride 
compound and the B^aTea film [11 S SI, SI . There- 
fore, the electron (hole) mobility fJ-^^l^i^-j along the trigo- 
nal axis is 300cm^/Vs (189cm^/Vs). From the classical 
definition of the carrier mobility and from Matthiessen's 
rule for the relaxation time, the electron (hole) mobility is 



expressed by /i 



[015] _ 



(cos^(0oMiI7+sin^(^')M7;;)- 



■ih) 



for the growth direction [015] 0, [s^]. It is equal to 
653cm'^/Vs (346cm^/Vs) for electrons (holes). In the 
next subsection, we present the temperature dependence 
of the bismuth telluride nanowire thermoelectric parame- 
ters calculated by means of the formulas considered here. 



B. Temperature dependence 

Fig. [5] shows the temperature dependence of the calcu- 
lated Seebeck coefficient for the intrinsic and the p— type 
bismuth telluride bulk material and the nanowires with 
different cross sections, growth directions, and excess hole 
concentrations. The temperature dependence of the in- 
trinsic nanowire Seebeck coefficient S is monotonic. For 
the p— type nanowire, the temperature dependence of S 
tends to be monotonic at a small nanowire cross section. 
In the intrinsic nanowire, the decrease of the cross sec- 
tion leads to an increase of the Seebeck coefficient mod- 
ulus. The intrinsic Seebeck coefficient is negative be- 
cause the Fermi energy tends to the conduction subbands 
owing to the difference between the hole and electron 
masses. For a high temperature, the calculated value 
of the intrinsic bulk Seebeck coefficient approaches the 
measured value -90 ^V/K, at 700 K [11]. The calcu- 
lated maximum value of the Seebeck coefficient S — 230 
/iV/K of the bulk p— type material shown in Fig. ^h) 
is less than the measured value S = 260 /xV/K at the 



room temperature [14|. This is because the Seebeck co- 
efficient was measured at the excess hole concentration 
Pex = 4 X 10^^ cm~^ which is less than the excess hole 
concentration p^x = 5 x 10^* cm~'^ used in our calcula- 



tion. In contrast to the p— type nanowire, the absolute 
value of the bipolar Seebeck coefficient for the growth 
direction [015] is less than that for the direction [110] 
owing to the carrier mobility difference as it follows from 
Eq. (29). In the expression (24) for the p— type See- 
beck coefficient, the hole mobility is reduced. An in- 
crease of the hole concentration leads to a decrease of 
the nanowire Seebeck coefficient. For example, S at 
Pex = 5 X 10^^ cm~^ is greater by 50 /xV/K than S at 
Pex = 1 X 10^^ cm""^ in the whole temperature range. In 
this case, the size quantum limit does not provide a large 
value of the p— type nanowire Seebeck coefficient. 

Excess holes suppress the difference between the values 
of the Seebeck coefficient corresponding to the nanowire 
growth directions [015] and [110]. In the tempera- 
ture range from 320 K to 420 K, S for the p— type 
nanowires with different cross sections and the growth 
directions are similar for the excess hole concentration 
= 5 X 10^^ cm~3. When temperature T > 300 K, the 
p— type bulk bismuth telluride material tends to be in- 
trinsic; as a result, the bulk Seebeck coefficient strongly 
decreases [ll] . Both the excess holes and the carrier con- 
finement shift the maximum value of the nanowire See- 
beck coefficient towards higher temperatures compared 
to that for the bulk material, because they suppress the 
intrinsic type of conductivity. For the excess hole con- 
centration Pex = 5 X 10^^ cm~'^, the maximum value of 
the Seebeck coefficient for the p— type square bismuth 
telluride nanowires with thickness 7 nm and the growth 
direction [110], as well as with thickness 15 nm corre- 
sponding to the growth directions [110] and [015], is equal 
to 250 fiV/K, 240 nV/K, and 230 nV/K, respectively. 
These values are greater than the respective bulk value 
of S, 230 /iV/K 13]. At the room temperature, the ab- 
solute value of the bipolar nanowire Seebeck coefficient 
is larger than 5* for the p— type nanowire with a small 
cross section. The maximum value of S for the intrinsic 
nanowire is achieved at low temperatures in contrast to 
that for the p— type nanowire. 

Fig.[n]represents the temperature dependence of the to- 
tal thermal conductivity k of the bismuth telluride com- 
pound as well as of the intrinsic and the p— type square 
nanowires. The thermal conductivity increases with tem- 
perature because the electron and hole parts of the to- 
tal thermal conductivity begin to play essential role at 
high temperatures. The bulk value of k measured along 
the cleavage plane is less than that for the 7 nm thick 
nanowire with the corresponding growth orientation at 
high temperatures. This difference is accounted for by a 
decreasing temperature dependence of the realistic car- 
rier mobility, which is supposed to be constant in our 
approximation. The maximum thermal conductivity of 
the p— type nanowire 0.42 W/Km (1.8 W/Km) is less 
than that of the intrinsic nanowire 2.5 W/Km (6 W/Km) 
with the same thickness 7 nm (15 nm). This increase is 
due to the carrier mass temperature dependence in the 
intrinsic nanowire. As we mentioned above, the phonon 
part of the nanowire thermal conductivity is less by one 
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FIG. 5: Temperature dependence of the Seebeck coefficient of (a) tlie intrinsic and (b) tlie p— type bismuth telluride compound 
(filled squares) and the square nanowire with thickness 7 nm (solid line), 15 nm for the growth direction [110] (dashed line) 



and [015] (dashed-dotted line) at pe 
direction [110] at Pex = 1 x 10^ 



5 X 10 cm . The dotted line corresponds to nanowire thickness 7 nm and the growth 
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FIG. 6: Temperature dependence of the thermal conductivity in (a) the bulk material (filled squares) fl^ j. the intrinsic square 
nanowire and (b) the p— type square nanowire with thickness 7 nm (solid line), 15 nm for the growth direction [110] (dashed 
line) and [015] (dashed-dotted line). The excess hole concentration is Pex = 5 x 10^* cm~^. For the growth direction [110], a 
dotted line corresponds to the nanowire with thickness 7 nm at pex = 1 x 10^^ cm~'^. 



order of magnitude than its bulk counterpart. There- 
fore, there is a great difference between the values of the 
thermal conductivity in the bulk and in the nanowire at 
low temperatures. In contrast to the Seebeck coefficient, 
the variation in the values of the thermal conductivity 
increases with temperature for the nanowires with differ- 
ent cross sections and growth orientations. The thermal 
conductivity of the nanowire with the growth direction 
[015] is less than that of the nanowire with the direc- 
tion [110] because of the greater effect of the confine- 
ment and another carrier mobility. The nanowires with 
smaller cross sections have a smaller thermal conductiv- 
ity. Hence, the carrier confinement causes a decrease of 
the nanowire thermal conductivity. For the cross section 
7x7 nm^, the thermal conductivity rises with increasing 
the excess hole concentration. 

Fig. [7] demonstrates the temperature dependence 
of the figure of merit of the bismuth telluride com- 
pound, the intrinsic and the p— type square nanowires 
with different growth directions, the cross sections 



7 X 7nm^ and 15 x 15 nm^, the excess hole concentra- 
tions Pex = 5 X 10^^ cm"'' and Pex = 1 x 10"^^ cm"''. For 
the intrinsic nanowire with the growth direction [110], the 
figure of merit achieves rather small values 0.1 and 0.12 
corresponding to the nanowire cross sections 15 x 15 nm^ 
and 7x7 nm'^ at temperature 500 K. They are less than 
the calculated value for the bulk, ZT = 0.17, because the 
bulk carrier concentration is greater than its nanowire 
counterpart. At temperature 260 K, the bulk figure of 
merit 0.03 is less than 0.07 associated with the 15 nm 
thick intrinsic nanowire due to the smaller nanowire ther- 
mal conductivity. When the temperature is below 360 K, 
ZT of the 7 nm-thick nanowire is smaller than that of the 
15 nm-thick nanowire, because the difference between the 
nanowire carrier concentrations is in favor of the latter 
case. Otherwise, when T > 360 K, this dependence in- 
verts: a higher value of ZT for the 7 nm-thick nanowire 
is provided by both a larger value of the Seebeck coef- 
ficient and a smaller value of the thermal conductivity. 
Therefore, in the middle of the temperature range under 
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FIG. 7: Temperature dependence of the figure of merit of (a) the bismuth telluride compound (filled squares), the intrinsic 
square nanowire and (b) the p— type square nanowire with thickness 7 nm (solid line) and 15 nm, with the [110] (dashed line) 
and [015] (dashed-dotted line) growth direction at pex = 5 x 10^* cm~^. Dotted line corresponds to nanowire with thickness 7 
nm grown along direction [110], at pex = 1 x 10^^ cm~^. 



consideration, the carrier concentration mainly influences 
the figure of merit of the intrinsic nanowire, while at high 
temperatures, the latter mainly depends on the Seebeck 
coefficient and the thermal conductivity. Since the carrier 
concentration and the Seebeck coefficient for the 15 nm- 
thick intrinsic nanowire grown along the [015] direction 
are less than those for the bismuth telluride compound 
and the nanowire with the other growth direction, its fig- 
ure of merit ZT is less than the ZT value of the bulk and 
the nanowire with the growth direction [110]. Hicks et 
al. 0, [13 also obtained such a relation from the calcula- 
tion for the bismuth telluride quantum wells and wires. 
At the excess hole concentration pex = 5 x 10^* cm~'^, 
the calculated bulk figure of merit 0.4 agrees with the 
value 0.5 reported by Hicks et al. at the room temper- 
ature. Otherwise, the figure of merit ZT — 2 for the 
7 nm-thick nanowire is considerably greater owing to 
the smaller thermal conductivity at temperature 300 K. 
The figure of merit for the p— type nanowire increases 
with decreasing the nanowire cross section. The See- 
beck coefficient and the thermal conductivity play the 
central role in the temperature dependence of the fig- 
ure of merit for the p— type nanowire in the temperature 
interval under consideration. For the growth direction 
[110], the maximum value of the figure of merit for the 
p— type nanowire is equal to 1.4, 1.6, and 2.8, correspond- 
ingly, at temperatures 310 K, 390 K, and 480 K and the 
cross sections 30 x 30 nm^, 15 x 15 nm^, and 7x7 nm^ 
{Pex = 5 x 10^^ cm^^). At the room temperature, the 
figure of merit equals 1.2, 1.3, and 1.7, respectively. The 
maximum value of ZT corresponds to the temperature at 
which the nanowire Seebeck coefficient has its maximum 
value, while the thermal conductivity begins to increase. 
This value increases and shifts towards high temperatures 
when the carrier confinement increases. Our calculations 
show that the optimal value of the excess hole concentra- 
tion is about Pex — 5 X 10^* cm^'^. Any deviation from 
this value leads to a decrease of the figure of merit of the 
p— type nanowire. 



Table IIIII contains the calculated thermoelectric pa- 
rameters and the available experimental data on the ther- 
moelectric properties of the p— and n— type bismuth 
telluride nanowires, superlattices, and bulk materials at 
the room temperature. The values of the Seebeck co- 
efficient and the thermal conductivity calculated for the 
p— type Bi2Te3 nanowire with thickness 7 nm agree with 
those measured on the p— type Bi^T / Sh2T superlat- 
tice with period 6 nm, while the nanowire electrical re- 
sistivity is twice greater than the superlattice resistivity. 
Since the effect of the one-dimensional confinement on 
the carrier concentration is less pronounced than the ef- 
fect of the two-dimensional confinement and our model 
does not take into account the surface scattering and 
phonon confinement effect, the calculated value of figure 
of merit disagrees with that obtained experimentally for 
the superlattice. It is difficult to make comparison with 
the experimental data on the nanowires because their di- 
ameter is greater than 80 nm when the confinement does 
not play an essential role. 



V. CONCLUSIONS 

The carrier confinement influences the thermoelectric 
properties of the bismuth telluride nanowires with thick- 
ness less than 30 nm. The confinement leads to an in- 
crease of the Seebeck coefficient as well as to a decrease of 
the thermal conductivity and the carrier concentration in 
both the intrinsic nanowires and the p— type nanowires. 
While the nanowire cross section decreases, the Femri 
energy, the carrier concentration, the thermal conductiv- 
ity, and the Seebeck coefficient change more significantly 
for the growth direction [015] than for the [110] direction 
because the carrier mass components are smaller for the 
former direction. In contrast to the Bi2Te3 nanowires, 
the 5*0.37^60.63 nanowires are intrinsic when the tem- 
perature ranges from 77 K to 500 K. Therefore, the de- 
pendence of the carrier effective mass on temperature 
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TABLE III: The calculated thermoelectric parameters for the p— type and the intrinsic bismuth telluride nanowires at the room 
temperature, in comparison with the experimental data on superlattices (SL), nanowires and bulk material. 





S 


K 


P 


ZT 


Dimensions/ 




ifiV/K) 


(W/mK) 






thickness /period 


n-type -Bio.54Teo.46 nanowire " 


-35 


1 


1.4 


0.02 


81 nm 


p— type Bio.46reo.54 nanowire" 


260 








100 nm 


n-type Bi2Te3/ Bi2Te2.83Seo.17 SL* 


-238 


0.94 


1.2 


1.4 


6.0 nm 


p-type BiaTeg/SfcaTes SL'"= 


165 


0.65 


0.53 


2.4 


6.0 nm 


n-type i3io.36Teo.64 bulk material** 


-150 


1.8 


1.2 


0.35 


15 mmxl5 mmx3 mm 


p— type BiaTea bulk material 


160 


2.7 


0.6 


0.22 


15 mmxl5 mmx3 mm 


Intrinsic Bio.37Teo.63 nanowire^ 


-153 


0.72 


11 


0.085 


15 nm 


p— type BiaTes nanowire*^ 


143 


0.53 


1.2 


0.94 


7 nm 



"Reference 3. 
'Reference 35. 

'^R. Venkatasubramanian, Phys. Rev. B 61, 3091 (2000). 

"^ReferenceflO. 
'This work. 



should be taken into account in the calculations of trans- 
port properties of such nanowires. Six equivalent valleys 
in the bulk bismuth telluride are split into two-fold and 
four-fold degenerate valleys for the nanowire growth di- 
rection [110] and into 3 doubly degenerate valleys for the 
direction [015]. Mainly, the excess holes and the car- 
rier confinement have opposite effects on the nanowire 
thermoelectric parameters. There is an exception, when 
both the excess holes and the carrier confinement sup- 
press the intrinsic type of electrical conductivity in the 
p— type nanowires at high temperatures owing to a weak- 
ening of the carrier concentration dependence on temper- 
ature. The size quantum limit can be used if the cross 
section of the p— type nanowire is less than 8x10 nm^ 
(6 X 7nm^ and 5 x 5nm^) at the excess hole concen- 
tration pex = 2 X 10^® cm~'^ {pex = 5 X 10^* cm~'^ and 
Pex = 1 X lO""^^ cm~^ correspondingly) in the restricted 
temperature range. The ZT value of the intrinsic 
nanowire does not exceed the maximal bulk value 0.16, 
until the thickness of the square nanowire is at least 
larger than 7 nm. For the growth direction [110], the 
maximum value of the figure of merit of the p— type 



nanowire is equal to 1.4, 1.6, and 2.8, correspondingly, 
at temperatures 310 K, 390 K, and 480 K and the 
cross sections 30 x 30 nm^, 15 x 15 nm^, and 7 x 7nm^ 
{Pex = 5 X 10^^ cm~^). At the room temperature, the 
figure of merit equals 1.2, 1.3, and 1.7, respectively. For 
the p— type Bi2Te3 nanowire, the figure of merit ZT 
reaches its maximum value when the excess hole concen- 
tration varies in the interval Pex = (4 -;- 8) x 10^* cm~'^. 
In general, the carrier confinement increases the figure 
of merit of the p— type bismuth telluride nanowires and 
shifts its maximum value towards high temperatures. 
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